data<-read.table("./AFFY/sample.miss",head=T,sep="\t")
h1<-hist(data$F_MISS,breaks=seq(0,1,0.01),plot=FALSE)
h1$breaks<-h1$breaks[2:101]
final<-cbind(h1$breaks,h1$counts)
data<-read.table("./CGPANEL/sample.miss",head=T,sep="\t")
h1<-hist(data$F_MISS,breaks=seq(0,1,0.01),plot=FALSE)
h1$breaks<-h1$breaks[2:101]
final<-cbind(final,h1$counts)
data<-read.table("./EXOME/sample.miss",head=T,sep="\t")
h1<-hist(data$F_MISS,breaks=seq(0,1,0.01),plot=FALSE)
h1$breaks<-h1$breaks[2:101]
final<-cbind(final,h1$counts)
data<-read.table("./IBC/sample.miss",head=T,sep="\t")
h1<-hist(data$F_MISS,breaks=seq(0,1,0.01),plot=FALSE)
h1$breaks<-h1$breaks[2:101]
final<-cbind(final,h1$counts)
colnames(final)<-c("missingrange","affy","cgpanel","exome","IBC")
write.table(final,"sample_missing.txt",sep="\t",col.names=T,row.names=F,quote=F)
